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Simplified  formulae  are  obtained  for  the  motion  of  solitons  in  nonuniform  media,  and  the  results  are 
illustrated  for  the  case  of  a  sinusoidally  varying  medium. 


Recently,  Chen  and  Liu1  presented  an  approximate 
theory  of  soliton  propagation  in  media  with  slowly  vary¬ 
ing  inhomogeneity,  within  the  context  of  the  adiabatic 
approximation.  They  considered  the  one-dimensional 
nonlinear  SchrtSdinger  equation  for  the  amplitude  q  in¬ 
cluding  an  inhomogeneity  term  On,  namely, 

(»  +  2lglJ-5»(*)^0(x><)  =  °  ,  (1) 


q(x,  t)  =  2t;  sech{2t)[x  +  4  51  +  m(l  >]} 

x  exp(-t{[2$f  t(l)Hx  +  tti(/))-4(f-if)i  +r(t)\). 

fo) 

If,  for  simplicity,  henceforth  we  consider  On's  chosen 
such  that  On(x0)  =  0,  then  we  may  identify 

-24  =  Vq  .  (6) 


where  l  is  time  and  x  is  position.  They  showed  that  to 
second  order  in  the  inhomogeneity,  the  adiabatic  ap¬ 
proximation  leads  to  the  solution 

q=u(y,  l  )exp[-iyp(l )  - i r(t  )J ,  y  =x  +  m(l ), 

P(f)  =  f  /(s)ds,  m(t)=  2  (  p(s)ds, 

*^0  •'0 

r(t)-  f  [e(,s)-p2(s)-f(s)Ms)]cls  , 


Inspection  of  Eq.  (5)  indicates  that  in  the  presence  of 
inhomogeneity  the  soliton  acquires  a  time-dependent  ef¬ 
fective  group  velocity  v,  frequency  u>,  and  propagation 
constant  k  of  the  form 

v-v0-m(t  )/t ;  k  =  k0~p(l)  , 
w  =  u>0  +  [r(t )  +  Ml  )p{l )  -  k0 m(t ) j  fi  , 

where  i>0= -45,  k„=-25,  and  u>0= -4(5'- - 1)2)  are  the 
values  corresponding  to  zero  inhomogeneity. 


/(s)s  ~  6n(x) 


X  =X  <*  ) 


e(s)=  &i(x(s))  — ,v(s)/(s), 


where  u  satisfies  the  nonlinear  Schrddinger  equation  in 
the  absence  of  inhomogeneity,  i.e., 

(i  It*  ly*  +2iwl2)“(>'^)  =  0  (3) 

and  .v,  the  classical  approximation  for  the  pulse  center, 
satisfies  the  equation 

r“’  fJx 

JI()  2 [Q-6n(x)il/2  "l  (4) 


At  first  glance,  it  appears  that  to  obtain  the  quantities 
in  Eq.  (7),  one  must  first  evaluate  the  integral  in  Eq, 
(4),  invert  the  result  to  obtain  x  =  x(t),  and  use  the  lat¬ 
ter  to  determine  /=/(!)  and  e  =  c(t ).  Finally,  one  would 
perform  the  series  of  repeated  integrations  over  t  in¬ 
volving/  and  e,  as  indicated  in  Eq.  (2).  In  this  note, 
we  indicate  how  explicit  expressions  may  be  obtained 
for  r,  m,  and  t  (and,  consequently,  v,  k,  and  to)  which 
do  not  involve  any  integrations  beyond  an  appropriate 
equivalent  of  Eq.  (4).  The  required  results  follow  by 
changing  the  integration  variable  in  Eq.  (2)  from  I  to  x ; 
a  transformation  commonly  utilized,  for  example,  when 
calculating  particle  orbits  in  electric  and  magnetic 
fields.  When  this  is  done,  one  finds  directly  that  (for 
6h(a„)  =  0] 


where  Q  is  an  appropriate  constant.  ( We  have  corrected 
certain  errors  in  Ref.  1,  with  respect  to  signs  in  the 
expressions  for  y  and  r(t)  in  Eqs.  (2),  and  a  factor  of  2 
in  the  integral  in  Eq.  (4)].  Combining  Eq.  (2)  with  the 
known  solutions  of  Eq.  (3),  one  obtains  the  single-soliton 
solution 


pU)~lv„-R(l).  m(t)-v0t  —x(t ), 

»•(/)=-»■(/)  + !?„/«(/),  (8a) 

.  ,  /*x  t<> 

«(<)*{ **-fri|x(/ ))}*'*.  win*  I  e»i(x)j,/2rf\. 

At 

(8b) 
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To  obtain  the  relation  between  rand/,  one  need  not  car¬ 
ry  out  the  integration  in  Eq.  (4),  but  may  instead  utilize 


™=t-t 

9  K 


(9) 


Since  W  is  a  functional  of  x(t ),  Eq.  (9)  may  also  be  in¬ 
verted  to  obtain  x=x(t ).  The  quantities  u,  ft,  and  u>  now 
take  the  simple  form 


v=x(t)/t,  k  =  R(l),  w  =  b)0  +  kv  -W(t)/t. 


(10) 


Note  especially  that  the  soliton  group  ,'elocity  v  involves 
no  further  integrations  beyond  either  the  one  in  Eq.  (4) 
or  (8b).  Moreover,  we  have  derived  the  very  simple 
and  perhaps  not  unexpected  result,  that  within  the  adia¬ 
batic  approximation,  the  soliton  envelope  becomes  a 
function  of  x  -7(/).  When  6 n~  0,  Y-*  v0t  and  the  envelope 
depends  on*  -u0t  as  expected. 

As  an  illustration  of  the  present  formulation,  we  con¬ 
sider  the  case  of  a  periodic  inhomogeneity  of  the  form 

6rt(*)  =  C  sin2ax  .  (11) 

A  variation  of  this  type  could  apply,  for  example,  to  the 
self -focusing  of  optical  pulses  in  media  with  a  sinus¬ 
oidally  varying  refractive  index.2  Then,  Eq.  (8b)  yields 

(*o  =  0) 

W(l )  =  a~lk0E  [a,  ax(t )] ;  <F=  C/fe|; 

(12) 

t  =~^=F[a,ax(t)](2ak0)'1  , 

where  F  and  E  are  incomplete  elliptic  integrals  of  the 
first  and  second  kind,3  respectively.  One  obtains  v  di¬ 
rectly  by  inverting  Eq.  (12)  to  yield  x  =x(t ),  which  in 
the  present  case  requires,  in  general,  numerical  com¬ 
putation;  q  and  w  take  the  form 

Accession  For 

~~NTIS  GliUI 
DOC  TAB 
Unannounced 
Justif  ic-tioij 


ft  =  ft0[l  -  a3  sin2ax(f )] 1  /2  , 
u>  =  w0  +  kv  -kaE(o,  ax)(at  )~l  . 


(13) 


For  the  special  case  a2=  1,  the  integrals  in  Eq.  (12)  are 
elementary1  and  one  is  able  to  express  a  explicitly  as 
a  function  of  t , 


x  =  at~1sin~iy{t),  y{t )  =  (e5  -  l)/(es  +1), 
s  =  2av0t , 


(14) 


whence 

v  =  sin‘‘y(t  )/<*!  >  *  =  2k0es/2/(es+ 1), 

(15) 

u)  =  tv0  +  kv  -k0y(t)/at  . 

Note  that  it  is  also  possible  to  obtain  closed  form  ex¬ 
pressions  for  v,  ft,  and  co  for  a2#  1  by  approximating 
E  and  F  with  an  appropriate  series  expansion.  How¬ 
ever,  it  is  clear  from  Eq.  (10)  that  one  must  proceed 
beyond  the  linear  approximation  for  x(x  ~  v0t )  to  obtain 
a  change  in  the  soliton  group  velocity. 

In  conclusion,  we  have  obtained  simplified  formulae 
for  the  effective  soliton  group  velocity,  frequency,  and 
propagation  constant  in  media  characterized  by  slowly 
varying  inhomogeneity,  and  have  illustrated  the  results 
for  the  c '**■'!  of  a  sin2ax  variation. 


'll.  H.  Chen  and  C.  S.  Liu,  Phys.  Fluids  21,  377  (1978). 

2See,  e.g.,  G.  B.Whitham,  Linear  and  Nonlinear  Waves  (Wiley, 
New  York,  197-1),  Sec.  16.4. 

3S.  M.  Selby,  editor,  CRC  Standard  Mathematical  Tables 
(Chemical  Rubber  Company,  Boca  Raton,  Fla.,  1968),  16th 
ed.,  p.  497. 

4L.  M.  Milr  Thomson,  in  Handbook  of  Mathematicl  Functions, 
edited  by  M.  Abramowitz  and  I.  Stegun  (Dover,  New  York, 
1965),  p.  594,  Eqs.  (17.4.21)  and  (17.4.25). 
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